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Abstract 

Using data collected with the BESII detector at e + e~ storage ring Beijing Elec- 
tron Positron Collider, the measurements of relative branching fractions for seven 
Cabibbo suppressed hadronic weak decays D° — ► K~K + , tt + tt^, K~ K + tt + tt~ and 
7T + 7T + 7r~7r~, D + — > K°K + , K~ K + n + and ir~n + n + are presented. 

PACS: 13.25.Ft, 14.40.Lb 



1 INTRODUCTION 



Hadronic decays of charmed mesons have been extensively studied. Measure- 
ments of relative lifetimes and semileptonic branching fractions for charmed 
mesons D + and D° suggest the presence of nonleptonic processes which en- 
hance the D° and suppress the D + width, and lead to the conclusion that the 
simple spectator model of charmed-meson decay is inadequate. As shown in 
Figure 1, all weak decays of heavy mesons may be described by six quark- 
diagrams: the external W-emission diagrams (a), the internal W-emission dia- 
grams (b), the W-exchange diagram (c), the W-annihilation diagram (d), the 
horizontal W-loop diagram (e), and the vertical W-loop diagram (f) [1]. Thus, 
a further understanding of the D decay mechanism, such as the contributions 
of other quark- diagrams and final states interactions, requires a systematic 
study of the hadronic decays. 
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(a) 



(c) 

c 

(e) 

Fig. 1. Typical Feynman diagrams for Cabibbo suppressed decays of charm mesons. 

At present, many experiments, such as, MARKII [2], MARKIII [3], E691 [5], 
E687 [6], E791 [4], FOCUS [8] and CLEO [7], have reported their measure- 
ments of Cabibbo suppressed hadronic decay fractions of D mesons. Our 
present measurements are based on a data sample of integrated luminos- 
ity of ~ 17-Spb- 1 at ^(3770) peak = 3.773 GeV) and ~ lG.Spb^ 1 for 
■^(3770) peak scan collected by Beijing Spectrometer (BESII) detector at e + e~ 
storage ring Beijing Electron Positron Collider (BEPC) [9]. This paper re- 
ports our measurements of Cabibbo suppressed relative branching fractions 
for several hadronic decay modes of charmed D mesons D° — > K~K + , 7r + 7r~, 
K~ K + n + n~ , 7T + 7r + 7r~7r~, D + — > K°K + , K~K + ir + and (through out 

this paper the charge conjugate states are implicitly included). 




2 BESII DETECTOR 



The Beijing Spectrometer (BESII) is a conventional cylindrical magnetic de- 
tector that is described in detail in Ref. [10]. A 12-layer Vertex Chamber (VC) 
surrounds the beryllium beam pipe and provides trigger information, as well 
as coordinate information. A forty-layer main drift chamber (MDC) located 
just outside the VC yields precise measurements of charged particle trajec- 
tories with a solid angle coverage over 85% of 4n; it also provides ionization 
energy loss (dE/dx) measurements which are used for particle identification. 
Momentum resolution of 1.7%yJ\ + p 2 (p in GeV/c) and dE/dx resolution for 
hadron tracks of ~ 8% are obtained. An array of 48 scintillation counters 
surrounding the MDC measures the time of flight (TOF) of charged particles 
with a resolution of about 200 ps for hadrons. Outside the TOF counters, a 
12 radiation length, lead-gas barrel shower counter (BSC), operating in lim- 
ited streamer mode, measures the energies of electrons and photons over 80% 
of the total solid angle with an energy resolution of a E /E = 0.22/v^ (E 
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in GeV). Outside the solenoidal coin, which provides a 0.4 T magnetic field 
over the tracking volume, is an iron flux return that is instrumented with three 
double-layer muon counters that identify muons with momentum greater than 
500 MeV/c. 

In this analysis, a GEANT3 based Monte Carlo package (SIMBES [11]) with 
detailed consideration of the detector performance (such as dead electronic 
channels) is used. The consistency between data and Monte Carlo has been 
carefully checked in many high purity physics channels, and the agreement is 
reasonable. 



3 EVENT SELECTION 



Charged tracks are required to satisfy \cos$\ < 0.8, where 9 is the polar angle 
in the MDC, and have good helix fit. The tracks that are not associated with 
Kg reconstruction are required to be originate from interaction point. Pions 
and kaons are identified by requiring the confidence level of desired hypothesis 
using combined measurements of time-of-flight [14] and energy loss in drift 
chamber to be greater than 0.1%. In addition, kaon and pion candidates are 
further identified by requiring the normalized weights, which is defined as 
CL a /(CL w + CL K ), where a denotes desired particle, exceeding 50%. 

K s candidates are detected through the decay of K s — > 7r + 7r~. Each oppositely 
charged track pair is assumed to be n + and tt~ . The decay vertex of Kg is 
required to be 5mm far away from the beam axis. The invariant mass 

is required to be within 0.020 GeV/c 2 of the Kg nominal mass. 

Cabibbo suppressed hadronic decay modes are expected at a lower rate (~ 
tan 2 8 C , where 9c is "Cabibbo angle") compared to relevant Cabibbo favoured 
modes, for which n/K misidentification becomes significant. The unique en- 
ergy of D meson at the ?/>(3770) can be exploited to reduce explicitly back- 
ground due to incorrect particle assignment. A single particle misidentification 
results in a reflection peak separated from the beam energy. 

The distributions of energy difference {AE) between measured energy of D 
candidates (E t3ig ) and beam energy (E^) are shown in Figure 2 for 4 Cabibbo 
allowed decay modes with correct n/K assignments and the reflection AE 
distribution with a single particle misidentification. AE of D candidate is 
required to be less than that of similar topological decay mode to reject particle 
misidentification combination. AE is further required to be within 50-100 MeV 
for different decay channels. 

The pair production of DD at ^(3770) provides a variable, which is defined 
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AE/GeV 



Fig. 2. AE distribution and reflection for similar Cabibbo suppressed modes by 
using the sample of Cabibbo favoured mode (a) D° — ► K~ir + (b) D° — > K°7r + (c) 
D+ — > K _ 7T + 7r + (d) D° — > K _ 7r + 7r + 7r _ . The distributions show the well-separated 
AE 1 peaks. 

as the beam-constrained mass 



4 = ^L-(Er) 2 

exploiting the fact that the total energy of all decay products must sum to 
the beam energy. As the uncertainty in the beam energy is much smaller than 
the uncertainty in the total reconstructed energy of the decay tracks, this 
approach yields much improved mass resolution compared to the invariant 
mass technique. 

The QED processes, t + t~ pair productions and cosmic backgrounds may con- 
tribute to D tags. Both of them have a lower multiplicity than that of DD 
decay, the requirement of N<± + N ncu /2 > 3 will eliminate most of these back- 
grounds, where N c ^ and N ncn represent the total number of charged tracks 
and neutral tracks respectively. At ^(3770), DD are produced with the angu- 
lar distribution sin 2 9 D) where 9d is the production angle of ^(3770) — > DD, 
|cos^d| < 0.8 is imposed to each D tag to enhance signal to background ratio. 

One event could be counted more than once as a tag candidate. In order 
to calculate the actual number of tagged events in an unbiased manner, the 
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following criterion is applied to select only one tag combination per event: if 
more than one combination of tracks form the desired tag, the combination is 
chosen when the lowest momentum track of this combination has the largest 
momentum of all other combinations. The result mass plots are shown in 
Figure 3, 4. 



4 DETECTOR ACCEPTANCE 



The detection efficiency is determined by a detailed Monte Carlo simulation 
of DD production, decay and detector response. The decay branching ratios 
of neutral and charged D meson are taken from the world average values [13], 
some unseen decay modes are set according to the rules of isospin conservation. 
Simulated events are processed through the event reconstruction, selection and 
analysis program. 

There are several sub-resonant decay modes in 3-body and 4-body Cabibbo 
suppressed channels. The detection efficiency is not uniform among these decay 
modes. The relative decay branching fractions and their errors in PDG are 
quoted. For D + — > it~it + -K + channel, p°7r + and 7r~7r + 7r + modes are considered; 
for D + — > K~K + tt + channel, 07r + , K*°K + and K~ K + n + are considered; for 
D° -> K-K+Ti+Ti- channel, ^Tr+vr", 0p°, K^K*°, K'K+p and R-K+ir+ir- 
are considered; for D° — > it~it + it + it~ channel, the uncertainty to Monte Carlo 
efficiency is estimated to be less than 2%. 



5 RESULTS 



The observed number of each Cabibbo suppressed decay channel is determined 
by fitting the distributions to a function of the form: 



F{m) = d\ 
+ 



( HL 



a 3 exp 



) exp \a 2 

(m-Mp) 2 
2a2 



where the first term parameterizes the background; the Gaussian term ac- 
counts for the signal. Just above the D mass, there is no more phase space 
available for a decay to a pair of DD mesons. The background term is AR- 
GUS form [12], where a\ is the normalization factor; a 2 , a scale factor for the 
exponential term; E^, the beam energy and fixed at 1.8865GeV while fitting 
the mass plot. For -^(3770) scan data, the beam energy is not a constant, but 
the background shape is similar, a constant background term 04 is applied to 
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Fig. 3. Beam energy constrain mass distributions for decays (a) D — ► K K + , (b) 
L>° 7T-7T+, (c) L>° -> and (d) D + -> (e) L>+ -> ^°vr + . 

evaluate the varying beam energy points. The mass resolution of each Cab- 
bibo suppressed mode is determined by fitting the Mt, c plot of similar Cabbibo 
allowed modes. The event number for each mode is summarized in Table 1. 



The pions that are part of reconstructed Kg's, are not used in modes ttty, 
7T7T7r, mnnT and KKw, KKinr. However, there remain some D decays into 
K s tc + , Kg7T + ir~ and K~K + K$, where the K$ is not identified as a separated 
vertex, which will be the major background source for decay mode D + — > 
7r~7r + 7r + and D° — ► 7r _ 7r + 7r + 7r _ , K~K + 7t + tt~ . To reduce the feed-down Kg — > 
7r + 7r _ background, a cut of \M n + n - — M K o\ > 0.040 GeV/c 2 is imposed on the 
invariant mass of each pion pairs for 7r~7r + 7r + and K~ K + it + tt~ modes. 

For the decay D° — > 7r _ 7r + 7r + 7r _ , plots of the invariant mass of all 7r + 7r~ 
combinations in D° candidates within the signal and sideband regions are 
shown in Figure 5(a) and (b) respectively. There is a clear Kg peak within D° 
signal region with a fitted events number of 112.5 ± 19.0 events and no clear 
Kg peak within the sideband region. The K s number is consistent with the 
expected background (Monte Carlo study gives this number as 98.8) due to 
D° — ► K°7v + 7v~ . These are thus subtracted from the in signal. 



6 SYSTEMATIC UNCERTAINTIES 



In this analysis, we normalize the relative ratios of Cabibbo suppressed decay 
to similar Cabibbo favoured modes, which permits the cancellation of many 
common systematic errors. The systematic uncertainties on energy cut and 
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Fig. 4. Beam energy constrain mass distribution for decay (a) D + — ► K~ K + ir + , 
(b) D+ -> 7T-7r+vr + , (c) D+ ,Fs:-7r + 7r + and (d) L>° -> K-R+it+it-, (e) 
£>° -» 7r-7r + 7r + 7r-, (f) L>° -> i^-vr+vr+vr-. 




C 1 1 1 3 

0.4 0.45 0.5 0.55 0.6 



Invariant Mass of 7C + it (GeV/c 2 ) 

Fig. 5. M n + n - for 7r _ 7r + 7r + 7r _ candidates within D° signal region (a) and sideband 
region (b). 

particle identification can not be canceled completely. 

Systematic uncertainties on particle identification and AE cuts are about 1-3 
% and 1-6% contributing to the relative fractions respectively. The systematic 
uncertainties on background subtraction due to the Kg contamination are 
estimated to be ~ 2% [15]. 

The detection efficiencies are not uniform among the different sub-resonant 
decay modes in the 3 and 4 body decays, 0.5 — 2.0% are estimated for different 
modes 7r~7r + 7r + , K~K + tt + and 7r~7r + 7r + 7r~. Large uncertainty in K~ K + n + n^ 
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mode is found due to the decays of low momentum kaons. 



7 CONCLUSIONS 



Table 1 

Measurement results of 7 Cabibbo suppressed decay modes. 



Decay mode 



Yield 



Relative efficiency 



Branching ratio 



K~K+ "~ 242.2±20.1 

K-TT+ 1934±49 

7T-7T+ 75.9±14.7 
K-TT+ 1934±49 

~K^K+ 63.2±9.8 
K°ir+ 287±18 

K-K+TT+ 181.2±20.2 

K-iT+ir+ 2324±53 

-K--K+-K+ 84.9±22.4 

K-TT+TT+ 2324±53 

K-K+ir+7r- 19.3±8.0 
K-ir+ir+ir- 1540±51 

TT-n+TT+n- (274.4±31.8)-(112.5±19.0) 

K-TT+TT+TT- 1540±51 



1.029 ±0.017 
1.146 ± 0.030 
0.991 ± 0.041 
0.669 ±0.015 ±0.010 
0.888 ±0.029 ±0.004 
0.286 ±0.021 ±0.017 
1.336 ±0.028 ±0.027 



0.122 ±0.011 ±0.004 
0.034 ± 0.007 ± 0.001 
0.222 ±0.037 ±0.013 
0.117 ±0.013 ±0.007 
0.041 ±0.011 ±0.003 
0.044 ± 0.018 ± 0.005 
0.079 ± 0.018 ± 0.005 



The relative fractions of seven Cabibbo suppressed decay modes are tabulated 
in Table. 1. For each mode, fitted event number, background number, rela- 
tive efficiency and relative branching ratio are listed. The errors of relative 
efficiency are Monte Carlo statistical error and systematic error due to sub- 
resonant modes respectively. The first error of branching ratio is statistical, the 
second one is systematic. The estimated systematic uncertainty is 3-6% for all 
modes, and the statistical uncertainties of the measurements are about 10% or 
greater. Results from this measurement are consistent with the world average 
values and we have improved the previous measurements of the D + — > K°K + 
relative branching ratio. The measurements of Cabibbo suppressed branching 
ratios presented here provide new insights into the mechanism of nonleptoinc 
D decays. Exact SU(3) symmetry predicts the equality of T(D° — > 7r~7r + ) 
and T(D° — > K~K + ). But the above results show they are not equal. Several 
distinct effects could contribute to this inequality. Final states interactions 
breaking SU(3) symmetry could however account for the difference. 
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